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Abstract 



A mesoscopic model of a diblock copolymer is used to study the stability of 
a lamellar structure under a uniform shear flow. We first obtain the nonlinear 
lamellar solutions under both steady and oscillatory shear flows. Regions of 
existence of these solutions are determined as a function of the parameters of 
the model and of the flow. Finally, we address the stability of the lamellar 
solution against long wavelength perturbations. 

I. INTRODUCTION 

We study a mesoscopic model of a block copolymer to describe the re-orientation of a 
lamellar structure by an imposed uniform shear flow that is either constant or periodic in 
time. This is a first step towards understanding known phenomenology pertaining to the 
response of the block copolymer microstructure to shear flows near the isotropic to lamellar 
transition The model that we use is based on the free energy of the diblock copolymer 
obtained by Leibler ||, and later by Ohta and Kawasaki ||, to which an advection term is 
added to incorporate the effect of the externally applied shear flow. We identify spatially 
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periodic solutions that correspond to a lamellar structure, and determine their stability 
against a number of long wavelength perturbations. 

Modulated phases are ubiquitous in physical and chemical systems [|J. They generally 
result from the competition between short and long range forces. Additional symmetries of 
the system (e.g., translational or rotational invariance) often lead in practice to rich textures, 
especially in systems of extent that is large compared with the characteristic wavelength of 
the modulation. Modulated phases often have interesting macroscopic behavior, and exhibit 
a complex response to externally applied forces. While it is possible to devise approximate 
constitutive laws to describe the macroscopic response of such phases, it is often neces- 
sary to explicitly address their evolution at the mesoscopic scale, and to determine how 
microstructure evolution influences the macroscopic response. 

We focus here on the lamellar phase observed in diblock copolymers below the order- 
disorder transition |||||5|]. Diblock copolymers are formed by two distinct sequences of 
monomers, A and B, that are mutually incompatible but chemically linked. At sufficiently 
low temperatures, species A and B would segregate to form macroscopic domains, but 
the chemical bonding between the two leads to a modulated phase instead. The detailed 
equilibrium microstructure depends on the relative molecular weight of the chains [f2],|3],|5],[7| 
and has been studied in detail within a mean field approximation fg,[| . 

We follow in this paper the approach of Leibler who introduced an order parameter 
field that describes the local number density difference of monomers A and B. The 
order parameter is defined to be zero above the order-disorder transition, and is finite and 
nonuniform below. Leibler's analysis was restricted to the weak-segregation limit (close to 
the order- disorder transition) within which the thickness of the interface separating the A- 
rich from the A-poor regions is of the order of the wavelength of the microstructure. Later, 
Ohta and Kawasaki extended Leibler's free energy to the strong segregation range, and 
showed the importance of long ranged effective interactions that arise from the connectivity 
of the polymer chains. We use this latter free energy as the driving force for the re-orientation 
dynamics, allowing also for passive advection of the order parameter by an imposed shear 
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flow. The model studied is similar to that considered by Fredrickson [Hj, except that we 



neglect thermal fluctuations and assume that both phases have the same viscosity. 

The stability of a lamellar structure to secondary instabilities has already been addressed 



in the literature, although in the absence of shear flow |lT| . |T2|j . In fact, the similarity be- 
tween the equations governing the motion of the lamellae and the Swift-Hohenberg model 
of Rayleigh-Benard convection [[i3H T5|| gives rise to a common phenomenology [12|. The 



lamellar structure is found to be stable only within a range of wavenumbers. At higher 
wavenumbers it undergoes an Eckhaus instability which generally results in a decrease of 
wavenumber, whereas for wavenumbers below that range the structure undergoes a zig-zag 
instability. In this paper, we extend these stability results to explicitly include fluid advec- 
tion by the imposed shear. We find that the stability boundaries are modified with respect 
to the zero velocity case in a way that depends not only on the amplitude of the shear, but 
also on the orientation of the lamellae relative to the flow. Of course, the latter dependence 
is absent in earlier treatments that neglected advection. 

Our results are a first step towards understanding the complex re-orientation phe- 
nomenology that has been observed experimentally [|TJ. In this initial analysis, we intro- 
duce a number of restrictive assumptions that we plan to relax in future work. First, our 
calculations are primarily two dimensional and thus can only address the so-called parallel 
and transverse orientations. Second, and more importantly, we neglect thermal fluctuations 
and any viscosity contrast between the two phases, elements that have been argued to be 
important in determining the main qualitative features of the re-orientation process. We 
also neglect flow induced by the lamellae themselves in response to the applied shear. These 
secondary flows could become important for the late stage coarsening of the lamellar struc- 
ture. Finally, we have confined our study to locating the boundaries of several secondary 
instabilities of the lamellar structure, but have not addressed the evolution following the in- 
stabilities, nor the coarsening of the resulting textured pattern p!6|JT7| . However, our results 
concerning the periodic base solution under flow, and its stability against long wavelength 
perturbations are the prerequisite building blocks of a more general theory. 



II. MESOSCOPIC MODEL EQUATIONS 



Following Leibler 0, we introduce an order parameter field, ip(r), function of the local 
density difference of monomers A and B. For a block copolymer with equal length sub-chains, 
the order parameter is ip(r) = PA<r ^ l ~ B ^ , where px,X = A, B is the density of monomer 
X, and po is the total density, assumed constant (incompressibility condition). A mean field 
free energy ^-"[^(r)] was derived by Leibler for a monodisperse diblock copolymer melt 0, 
and later by Ohta and Kawasaki 0. In units of fceT, where ks is Boltzmann's constant and 
T the temperature, the free energy is comprised of two terms, T j (p Q k B T) = T s + T\. The 
term T s incorporates local monomer interactions, 

?s = /drg|V^| 2 -^ 2 + ^ 
and is formally identical to the Ginzburg-Landau free energy commonly used to describe 



phase separation in a binary fluid mixture [fLq| . The parameters k, t and B can be ap 



proximately related to the polymerization index N, Kuhn's statistical length b and the 



Flory-Huggins parameter \ through the relations n = T = 2xN m ~ an d B = 4|t2 111 



3 ' n a 1KX N 2 b 2 

Long range interactions arising from the covalent bond connecting the two sub-chains 
are contained in T\, 

f l = ?-JJ dr dr'G(r - r')^(r)^(r') 

where the kernel G(r — r') is the infinite space Green's function of the Laplacian operator 
V 2 G(r — r') = —5(r — r'). The nonlocal interactions arising from the connectivity of the 
chains lead to a thermodynamic equilibrium state with a nonuniform density. In our case of 
equal length sub-chains, the equilibrium configuration is a periodic lamellar structure, with 
a characteristic wavelength of the order of 100 A for a typical system. 

Given this mean field free energy, a phenomenological set of equations that govern the 
temporal relaxation of equilibrium thermal fluctuations of ip{ r ) an d of fluid velocity v has 



been derived close to the order- disorder transition [ 19, 20, 10 1. A similar phenomenological 
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description can be used below the order- disorder transition under the assumption that the 
local relaxation of the order parameter field at the mesoscopic scale is still driven by mini- 



mization of the same free energy f21|]2!q| . Under this assumption, ip obeys the time-dependent 
Ginzburg-Landau equation, 

* + v.V, = *vf, (!) 

where M is a phenomenological mobility coefficient, and S /Sip stands for functional differ- 
entiation with respect to ip. Equation flU) includes the effect of advection by a local velocity 
field v, which satisfies an extended Navier-Stokes equation 

+ v • V v = z/V 2 v - + 2 

at p dip p 

where v is the kinematic viscosity of the fluid, assumed constant and independent of ip, 
p is the fluid pressure, and appropriate boundary conditions for both ip and v must be 
introduced. The last term on the right-hand side of Eq. (fj) is required to ensure that there 
cannot be free energy reduction by pure advection of ip . This term is sometimes referred 
to as osmotic stress, and it leads to the creation of rotational flow by curved lamellae that 
is directed towards their local center of curvature. 

We focus on a layer of block copolymer, unbounded in the x and y directions, and 
being uniformly sheared along the z direction (Fig. [I]). The layer is confined between the 
stationary z = plane, and the plane z = d which is uniformly displaced parallel to itself 
with a velocity f p i an e = sd in the case of a steady shear, and f p i an e — jduj cos(ut) in the 
case of an oscillatory shear, s is the dimensional shear rate in the steady case, and 7 is the 
dimensionless strain amplitude in the case of an oscillatory shear of angular frequency uj. 

The general problem defined by Eqs. ([I]) and (0) can be considerably simplified by noting 
that under typical experimental conditions inertia is negligible (u>d 2 jv ^ 1). Furthermore, 
we will neglect in this paper the term (5J r /5ip)'Vip/p in Eq. (|2|). Under these conditions, 
Eq. (Q) admits a simple solution that satisfies the specified conditions at the moving plane: 
v = s z i for a steady shear, and v = 7 d 00 cos(ut) z i for an oscillatory shear, where i is 



the unit vector in the x direction. Therefore, the problem reduces to a single governing 
equation for the order parameter field ip (Eq. (|I|)) under a prescribed advection velocity v. 
As discussed in the introduction, previous theoretical work on the formation and stability 
of lamellar structures further neglected advection of ip in Eq. ([[]). The results presented in 
this paper are free of this restriction. 

Since the base state to be considered is comprised of spatially uniform lamellae advected 
by the shear flow, it is convenient to introduce a new frame of reference in which the 
velocity vanishes. Define a new system of non mutually orthogonal coordinates (xi,X2,xs) 
by x\ = x — a(t)z,X2 = y and x% = z. The dimensionless quantity a(t) = st for a steady 
shear, and a(t) = 7sin(a;t) for an oscillatory shear. All the calculations reported in this 
paper, both analytical and numerical, have been performed in this new frame of reference. 
Analytical calculations consider an unbounded geometry in the X\ and £2 directions and 
periodic boundary conditions along the X3 direction, whereas the numerical computations 
have been conducted in a two dimensional, square domain on the (xi, £3) plane and consider 
periodic boundary conditions along both x\ and X3. Note that both frame of references 
coincide at t = 0, and at equal successive intervals of one half the period of the shear in the 
case of oscillatory shear. 



Dimensionless variables are introduced by defining a scale of length by J/c/r, a scale 
of time by k/Mt 2 , and an order parameter scale by ^Jt/u. In the transformed frame of 
reference and in dimensionless variables, Eq. (|l|) reads, 

^ = V' 2 (^ + ^ 3 -V^)-^, (3) 

with 

d 2 „ , , d 2 d 2 d 2 



7/2 



l + a\t) - 2a(t) _ _ + ^ + 



dx 2 dxidx s dx 2 dx\ 

There is only one dimensionless group remaining Bk/t 2 , which will be simply denoted by 
B in what follows. 



We will first show in Sec. [ITT] that below (but close to) the order-disorder transition point 



(in the weak-segregation limit), Eq. (|3|) admits periodic solutions. Their stability against 
infinitesimal long wavelength perturbations is the subject of Sec. fVj. 

III. LAMELLAR SOLUTION IN THE WEAK-SEGREGATION LIMIT 

In the absence of shear (a(t) = 0) the uniform solution of Eq. (Q), t/j — 0, loses stability 
at the order-disorder transition. In a mean field approximation the transition occurs at 
B c = 1/4. This is a supercritical bifurcation with a critical wavenumber q c = ^/l/2. Near 
threshold (0 < e = (B c — B)/2B C <C 1) there exist periodic stationary solutions of the form 

ip(r) = 2A cos(q • r) + A 1 cos(3q ■ r) + . . . , (4) 

with A 2 = q - ~ 0(e), and A\ of higher order in e. This solution only exists for a range 
of wavenumbers q such that u(q 2 ) = q 2 — q A — B > 0. 

For nonzero shear, we seek solutions of Eq. (|3|) of the form of Eq. (^), with r = (x±, X2, x^) 
expressed in the sheared frame basis set {ei = i, e 2 = j, e 3 = a(t)i + k}. Wavevectors are 
expressed in the reciprocal basis set {gi = i — a(t)k,g2 = j, g3 = k}. Therefore we keep 
the same functional form as for nonzero shear, but allow a time- dependent amplitude A(t). 
Note that the components of the wavevector q are assumed to be independent of time and 
given by q\ = q x {t — 0), q% — q y (t — 0) and g 3 — q z (t — 0) respectively. The wavevector itself 
depends on time through the time dependence of the reciprocal basis set. Such a solution 
corresponds to a spatially uniform lamellar structure with a time- dependent wavevector that 
adiabatically follows the imposed shear in the laboratory frame (see Fig. [1]). Inserting Eq. 
(|) into Eq. (|j), we find to order e 3//2 (a is itself of order e), 

A-3q 2 (t)A 3 , (5) 

with q 2 (t) = q\ + (a(t)gi — g 3 ) 2 + g|, and cr(q 2 ) = q 2 — q A — B. This nonlinear equation with 
time-dependent coefficients can be solved exactly in the two cases of steady and oscillatory 
shear flow. 



a 



7 



In the case of a steady shear a(t) = st. We find, 

( e 2H(t) w 1 -V2 



A(t) = + 6e 2 ^ / o ^e- 2 ^')g 2 (t')] , (6) 



with 



H(t) = (4 + B- q i)t + {i~ 2 q i)s qiq , e + Mg!M±i) t 3 _ g?g3S 3 1 4 + 5j£ (7) 



The constant quantity q = J q\ + q\ + g| is the initial wavenumber, and A(0) is the initial 
amplitude. For the special case q\ = 0, A(t) simply relaxes to its equilibrium value in the 
absence of shear A = " M — . This corresponds to an initial orientation of the structure 
which has no component transverse to the flow. For any other initial orientation, the shear 
induces changes in the lamellar spacing in the laboratory frame of reference (Fig. 0). As a 
result, the amplitude A(t) decreases and approaches zero at long times. Hence, the structure 
melts and reforms with a different orientation which we cannot predict on the basis of our 
single mode analysis. The emerging structure presumably results from the amplification of 
thermal fluctuations near the point at which the amplitude A(t) vanishes, and they have 
been neglected in our treatment. Thermal fluctuation effects have been accounted for by 
others 

For an oscillatory shear a(t) = jsm(ut). We first examine the stability of the uniform 
solution ip = against small perturbations. Linearization of Eq. (|5|) leads to, 

dA(t) 



dt 



a 



q\t) A(t), (8) 



with a(t + T) = cr(t) and T = 2-k/uj. Equation (H) constitutes a one-dimensional Floquet 
problem. The solution A = is unstable when 

a = f T a(t)dt > 0. (9) 

J 

The resulting neutral stability curve is given by, 

B - % ~ % g 2 ' ^ ' 



Instability modes can be conveniently classified by considering the relative orientation of 
the lamellae at t — and the shear direction. We define a parallel orientation, 93 7^ 0, 91 = 
92 = 0; a perpendicular orientation, q 2 7^ 0, 91 = 93 = 0; and a transverse orientation, 
qi 7^ 0, q2 = q3 = 0. The following instability points are identified depending on the 
orientation of the critical wavevector: a transverse mode with 



b - 1 ( 2 + ^ 2 ) 2 a - I 4 + 2 ^ 2 rm 

c 28 + 8 7 2 + 3 7 4 ' gic V 8 + 87 2 + 3 7 4 ' 1 j 



a mixed parallel-perpendicular mode with 
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B c = ^, q u = 0, 2q 2 2c + 2ql=l, (12) 



and a mixed parallel-transverse mode defined by 



1 7 7 2 + 16 / 1 fl / 37 2 + 8 

4 15 7 2 + 16' q2c ' gic V IV + 16' g3c V 30 7 2 + 32' 



Note that the threshold corresponding to perturbations of wavevectors that do not have a 
projection along the transverse direction are not affected by the shear. Furthermore, neither 
the stability boundaries nor the values of the critical wavenumbers depend on the angular 
frequency uj. 

In what follows, we consider mainly two dimensional solutions in the plane q2 = 
(transverse and parallel orientations) to make contact with two dimensional numerical cal- 
culations. As an example, Fig. ^| shows the neutral stability curve in the (91,93) plane for 
mixed parallel-transverse modes at e = 0.04, and for several values of the dimensionless 
strain amplitude 7 . Recall that 91 = q x (t = 0) and 93 = q z (t = 0) define the initial orienta- 
tion of the lamellae. The figure shows that the shear does not modify the neutral stability 
curve in the vicinity of 91 = (parallel orientation), whereas the curve is shifted near 93 = 
(transverse orientation). Large changes are observed for oblique wavevectors, including the 
complete suppression of the instability at sufficiently large values of the strain amplitude. 

Above threshold, Eq. (|5p can be solved to yield the time-dependent amplitude of the 
lamellar structure under oscillatory shear. We find, 



A(t) 



3 2(/(t)-c 4 -c 5 ) 



A(oy 



+ 6e 2/ W f dt'e- 2I ^q 2 (t') 
Jo 



-1/2 



The function I(t) is given by 



I(t) = C\t + c 2 sin(2c<jt) + c 3 sin(4c<jt) + c 4 cos(W) + c 5 cos 3 (u;£), 



with 



Cl 



3^ + (ggg±MzW + + b - g 2 



(14) 



(15) 



c 2 = - 



9 h 4 +(2gg+4gg-l)7 2 gf 
Au) 



C 3 



4 4 
32w 



C4 



( 4 gp -2)791 g3+47 3 gf 93 



and c 5 = -^ 



We note that the stability condition Eq. fllPf) is equivalent to c\ = 0. Hence, the 
asymptotic behavior of A(t) at long times changes qualitatively depending on the sign of 
c\. For Ci > 0, lim^oo e~ 2/ = 0, so that the integral in Eq. fll~4] ) tends to a finite constant. 
Since the prefactor e 2/ diverges exponentially, A(t) decays to zero. If, on the other hand, 
Ci < 0, A(t) becomes periodic at long times. To prove this statement, we first rewrite the 
second term inside the curly brackets in Eq. fll~4] ) as 

(16) 



T{t) =6e 2/ W f dt'e- 2c ^-^- 2 f {t '\ 2 (t') 
Jo 



with f[t') = c 2 sin^co't') + c 3 sin(4ut') + c 4 cos(u;t / ) + c 5 cos 3 (u;t / ). Since both f(t') and q 2 {t') 
are periodic with period T = 2tt/uj, we can decompose Eq. ( |T6"D into 



l(t) = 6e 2/ « 



V e 2jclT / dt' e -^'-m'+t )q 2 {t , + t) + e 2 Cl t 

3=1 J ° 



t-nT 



dfl e -2c lt '-2f(t') q 2(t') 



(17) 



where n is an integer such that < t — nT < T. In the limit of large t and with ci negative, 
the last term on the right-hand side of Eq. (|T7| ) vanishes while the sum X)?=i e 2 ^ ClT converges 
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to l/(e 2ciT — 1). Combining Eqs. ( ^4|) and (|P7D yields an asymptotically periodic solution 
for A(t), 

-1/2 

(18) 

The condition c\ = can also be understood in terms of a critical strain amplitude 7 C 
above which an existing lamellar structure of a given orientation at t = will melt (i.e., 
A(t) will decay to zero at long times) The value of 7 C that corresponds to c\ = is given by, 

7c = [(_ 6 + y/W-Adc)^] 1 ' 2 , (19) 

with b = {2ql + 4g| — l)q 2 /2, c = q^ + B — q^ and d = 3g^/8. Note again that the critical 
strain amplitude is independent of the angular frequency to. 

In order to test the approximations involved in Eq. (|]), namely that the wavevector q 
adiabatically follows the flow, and the single mode truncation for small e, we have undertaken 
a numerical solution of the model equation in a two dimensional, square geometry (see the 
Appendix for the details of the numerical method). As a first example, we consider an 
oscillatory shear of angular frequency ui = 0.02 imposed on a lamellar structure of initial 
wavevector (91,93) = (0.687,0.098). The critical strain amplitude for this initial orientation 
is 7 C = 0.695. Figure ||] shows the temporal evolution of A(t) for two values of 7, one larger 
and one smaller than j c . The solid lines are the predictions of Eq. flll]) , and the symbols 
are the results of the numerical calculation. The agreement in both cases is excellent. 

IV. SECONDARY INSTABILITIES OF THE LAMELLAR PATTERN 

In order to address the stability of the lamellar pattern, we next consider long wavelength 
perturbations of the base state with wavevector Q = (Qi, Q 2 , Q3) such that its components 
are also constant in the sheared frame of reference. Close to threshold, perturbations evolve 
in a slow time scale compared to the inverse frequency of the shear. We therefore assume 
that the wavenumber of any long wave perturbation would adiabatically follow the imposed 
flow. Specifically, we consider a solution of the form, 
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A(t) 



-2 Cl T _ I 



dt>e -2c lt '-2f(t' + t) q 2^ + ^ 



V>(r, t) = [A{t) + SA + e lQ - r + 5A_e^ Q - r ]e iq - r + c.c. 



(20) 



where A(t) is the nonlinear solution obtained in Section |TTl[ Substituting Eq. fl2~0"[) into Eq. 
(H), and linearizing with respect to the amplitudes 5A + and 6A_, we find, 



d_ 

Of 



5A + 




5A + 




= L(t) 




5A_ 




5A_ 



(21) 



with 



and U 



L(t) 



-(q±Q) : 



l 2 + -B + QA(t) 2 L 
3A(t) 2 L 



3A(t) 2 l + 
-l_-l 2 _-B + 6A(t) 2 L 



'l+a{t) 2 ){q l ±Q l ) 2 + 2a{t){q 1 ±Q 1 ){q 3 ±Q 3 )~{q 2 ±Q 2 ) 2 -{q 3 ±Q 3 



In general, the matrix elements are complicated functions of time, and we have not 
attempted to solve Eq. (pi]) analytically. For 7 < 7 C the operator L contains terms that are 
both periodic in time and decaying transients. At long enough times, A(t) is given by Eq. 
(fjjl), and the linear system Eq. (|2~ID has periodic coefficients. Hence, it reduces to a two 
dimensional Floquet problem for the amplitudes SA + and 5A^ [24| . 



In order to gain some insight into the stability problem, we first briefly review the known 



results for zero shear fTT| , |T2|| . In this case A(t) is a constant, and the matrix elements of L(t) 
are independent of time. An eigenvalue problem results by considering solutions of Eq. ( PT| ) 
of the form 5A± ~ e a+t + e cr - t , and instability follows when either eigenvalue is positive. Two 
modes of instability are obtained: a zig-zag (ZZ) mode that leads to a transverse modulation 
of the lamellae (Q ■ q = 0), and an Eckhaus (E) mode that is purely longitudinal in nature 
Q • q = Qq. In the zig-zag case, cr+(Q) has a maximum at 

1 - 2q 2 - 3A 2 



Qmax,ZZ 



(22) 



The eigenvalue cr + (Q) changes sign on the line q = q c , which therefore defines the zigzag 
stability boundary. In the Eckhaus case, we find after some straightforward algebra that the 
perturbation with the largest growth rate is 
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Qmax.E 



QA5q 4 



A5q 



2\2 



(23) 



64 5q 2 

with 5q = q — q c . Therefore the Eckhaus stability boundary is given by e = 12 5q 2 . These 
results are schematically summarized in Fig. |4j. The hatched area is the region of stability 
of a lamellar solution in the absence of shear flow. It is worth pointing out that this stability 
diagram is identical to that of the Swift-Hohenberg model of Rayleigh-Benard convection 
T5[ . Shiwa [|12| has recently shown that in the weak- segregation limit (e C 1), and in the 



absence of shear flow, the amplitude equation describing slow modulations of a lamellar 
solution is the same as the amplitude equation of the Swift-Hohenberg model near onset 
of convection. The same stability diagram has been derived by Kodama and Doi [JTTJ] by 
examining free energy changes upon distortion of a lamellar pattern. 

We now return to the Floquet problem of Eq. ( |2T| ) when A(t) is a periodic function of 
time (Eq. (|T8|)). Since A(t + T) = A(t) (T = 2n/u), the solution of (|2T|) is given by 



5A + 






5A_ 







with <f)±(t + T) 
(0,T), 



(24) 



(f)±(t). Equation fl2T|) is then transformed to an eigenvalue problem within 



d 








+ L(t) 




dt 


= —a 






4>4t) 







(25) 



Given that the function A(t) is quite complicated, we have solved this eigenvalue problem 
numerically. The eigenvalue a can depend in principle on the wavevector of the base state 
q, on the wavevector of the perturbation Q, and on the amplitude 7 and frequency u of the 
shear. For ease of presentation, we have focused on the case e = 0.04 although extension to 
other values of e is straightforward. 

Figures [5] and § summarize our results for the cases 7 = 0.2 and 7 = 0.4 respectively, 
and show the stability boundaries in the plane (qi, q^), as well as the neutral stability curve 
already shown in Fig. |2|. As before, (91,93) is the wavevector of the lamellar structure at 
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t = 0. At fixed e, 7 and cu, these curves have been obtained by determining the loci of q at 
which the function cr(Q) changes from a maximum to a saddle point at Q = 0. First we note 
that any orientation of the lamellar pattern that is not initially close to either parallel or 
transverse is unstable to moderate shears. Second, and up to our numerical accuracy, these 
curves appear to be independent of angular frequency. Finally, and contrary to the case 
of no shear, the reciprocal basis vectors are not time independent. Since the components 
of both q and Q are independent of time in the sheared frame, their mutual angle is not 
(except for the case in which they are parallel). However, the following statements can be 
made about the type of secondary instability. We have found that the secondary instability 
is of the longitudinal type only when either qi = or q 3 = (intersections between the lines 
marked with circles and the axes in Figs. |5] and ^|). Otherwise, the angle between q and Q 
is time- dependent. The lines on both figures marked with squares have the property that 
even though both q and Q are functions of time, their angle oscillates periodically around 
90°. 

The cases discussed up to now concern long wavelength instabilities of the base periodic 
pattern that are associated with the broken translational symmetry of the original system by 
the appearance of a periodic pattern. We now show that it is possible to obtain analytical 
expressions for the stability boundaries against finite wavelength perturbations that may 
have some experimental relevance as well. In some experimental protocols, the lamellar pat- 
tern is first obtained in the absence of shear. The resulting configuration comprises regions 
or domains of locally parallel lamellae but with a continuous distribution of orientations. 
A shear flow is then initiated and the reorientation of the pattern studied as a function of 
time. The pattern obtained in the absence of shear may be now unstable to several finite 
wavenumber perturbations that would not have been observable in the case in which flow is 
present throughout the ordering process. In the latter case the unstable orientations would 
have decayed away during the process of formation of the lamellae. In addition, the approx- 
imation that we derive below is generally valid when Q 3 cannot approach zero, as is the case 
in a system of finite extent in the direction of the velocity gradient. 
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We first define the following linear transformation, 







<5_ 





3A(t) 2 l- a + -a.+L 2 2-L 11 
3A(t) 2 l- a + -a-+L 11 -L 2 2 



a+—a— 2(<x+— cr_) 

which diagonalizes matrix L(t), and where 





8A+ 




5A + 






= T(t) 






<L4_ 




5A„ 



(26) 



<7±(f) 



L n + L 22 ± J(L n - L 22 ) 2 + 3Ql + l.A(t) 4 



(27) 



Combining Eqs. (|2l"l) and (p6j ), we find 
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ct+ 




5+ 




*-(*) 




ct_ 







yl(t) 2 L 



OA 



<7_ 



Mi -M 2 
-Mi M 2 



5_ 



(28) 



where Mi 



d_ 

M 



2A 2 l- 



and M 2 



9i 



2A 2 l. 



For finite Q, « L 

Assuming L n — L 22 > (the other case leads to no extra complications), a+ = L u and 
<t_ = L 22 . Also Mi = Ll ^i~^ 22 and M 2 = so that the equation for 5 + decouples from the 
equation for 5_. The solution for 5 + is 



5+(t) = 5 + (0)e J 



at' 



In 



L22— Ln 
A 2 l- 



(29) 



)df = ffc+d? = 0. We 



The stability boundary is defined by a = f (a + — ^ hi 
have checked that this stability condition agrees with the numerical stability analysis based 
on Eq. Q2§) for finite Q. 

We finish by illustrating the re-orientation dynamics of the lamellar structure following a 
long wavelength instability by direct numerical solution of the governing equation. We focus 
on the region in which the uniform lamellar structure is linearly unstable. The first example 
discussed concerns a lamellar structure of initial wavenumber (91,93) = (—0.4908,0.4908) 
being sheared periodically with an amplitude 7 = 1. Figures 0A,B,C show the sequence of 
configurations obtained when u = 5 x 10~ 6 . A long wavelength transverse modulation of 
the lamellae is observed (Figure |7]A). Subsequent growth leads to the formation of a forward 
kink band similar to that recently observed experimentally (Figure |7|B) p5,26|. As the 
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strain grows larger, the kink band disappears leaving behind a lamellar structure without 
any defects and oriented differently relative to the shear (Figure 0C). 

In the second example (Figs. |7p,E,F), a structure initially transverse to the flow is being 
sheared periodically with an amplitude 7 = 1 and a frequency uj = 5 x 10~ 7 . A longitudinal 
perturbation is clearly visible that manifests itself by local compression and dilation of the 
structure, leading to the disappearance of a pair of lamellae. The overall result is an increase 
in the lamellar spacing without any change in the orientation. 

In summary, we have obtained a nonlinear solution of the model equations that govern 
the formation of a lamellar structure in the weak-segregation limit. The solution is a pe- 
riodic lamellar structure with a time-dependent wavevector that adiabatically follows the 
imposed shear flow, and a time-dependent amplitude which we have computed for the cases 
of steady and oscillatory shears. In the case of an oscillatory shear, the periodic solution only 
exists for a range of orientations of the lamellae relative to the shear direction. The width 
of the region depends on the shear amplitude but not on its frequency. Long wavelength 
secondary instabilities further reduce the range of existence of stable lamellar solutions. The 
corresponding stability boundaries depend again on the shear amplitude, but are indepen- 
dent (up to our numerical accuracy) of frequency. We next plan to examine the stability of 
the nonlinear solution presented in this paper when neither osmotic stresses nor viscosity 
contrast are neglected. 
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APPENDIX A: NUMERICAL ALGORITHM 



We use a pseudo-spectral technique to solve Eq. (|) in two spatial dimensions and in the 
sheared frame with periodic boundary conditions along both directions. Equation ([3]) can 
be written as, 



^ = a(t)^-q 2 (t)^, 



(Al) 



where a{t) = q 2 {t) — g 4 (t) — B and q 2 (t) = q\ + [a{t)qi — q 3 ] 2 . The algorithm we use, due 
to Cross et al. J27|, is obtained by first multiplying both sides of Eq. flA~l]) by exp(— cr(t') t') 
and integrating over t'. This gives 

rt+At 



exp(-cx(t) t')i/j 



\t+At_ 



dt'^ 3 {t')exp{-a{t)t'), 



(A2) 



where we have assumed er(t') ~ cr(t) and q 2 (t') ~ q 2 (t). Next, we write the non-linear term 
i/j 3 (t') as a linear function of t' in the interval t < t' < t + At, i.e., 



»(Q W ^(t) + ^ + A ^-^(t-Q. 



At 



(A3) 



Combining the last two equations finally yields 

^(t + At) = exp(a(t) At)^(t) - q 2 (t)ip 3 (t) 
V 3 (t + At) - ^~ 3 (t) 



At 



exp(a(t) At) - 1 
exp(>(t) At) - (1 + cr(t)At) 



(A4) 



Eq. (|A4[) is first evaluated with the last term on its right-hand side set to zero. The resulting 
value for ip{t + At) is then used to estimate ip 3 (t + At). Eq. (|A4|) is finally applied a second 
time with all three terms on its right-hand side now included in the calculation. The fact 
that the nonlinear terms are integrated using an explicit procedure in time limits the size of 
the time step At that can be used in simulations of the model. 

All the numerical results presented in this paper were obtained in the sheared frame of 
reference with 128 x 128 spectral modes. We have chosen B = 0.23 (which corresponds to 
e = 0.04) and a time step of maximum size At = 0.2, for which no numerical instability was 
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observed. The initial condition ^(r, t = 0) is, unless otherwise noted, a lamellar structure 
obtained by numerical integration of Eq. ( |A4| ) with a(t) = (no shear) starting from a 
random initial condition (a gaussian distribution for ip if zero mean and small variance), for 
approximately 300,000 iterations until a stationary lamellar structure is reached. 
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FIGURES 



------ V 




sdt 



FIG. 1. Schematic representation of the configuration studied. We also show a schematic of the 
distortion of a lamellar pattern under uniform shear flow. The imposed velocity field is along the 
x direction. The velocity is specified at the z = d boundary, and vanishes at z = 0. The lamellae 
in this graph are transverse to the flow at t = (solid lines). At a later time (dotted lines) the 
lamellae are at an angle with respect to the flow, and the wavelength has changed accordingly. 
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FIG. 2. Neutral stability curves at e = 0.04 as a function of the wavenumber of the perturbation 
at t = (91,^3) for the values of the dimensionless strain rate 7 indicated. The inner regions 
bounded by the various curves represent the regions in which the solution ip = is linearly unstable. 
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FIG. 3. Temporal evolution of the amplitude A(t) of the single mode solution given in Eq. 



(14) for an oscillatory shear (solid lines), along with the corresponding numerical solution of the 
full model for e = 0.04, u = 0.02 and 7 = 0.5 (•), and 7 = 0.75 (o). Time has been scaled by 
the period of the shear T = 2tt/uj. For 7 = 0.75, the initial periodic pattern is unstable against 
uniform melting and the amplitude A(t) decreases to zero. On the other hand, for 7 = 0.5 a 
spatially periodic solution is stable, and after a transient, A(t) becomes a periodic function of time 
(Eq. ©). 
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q-q c 



FIG. 4. Schematic representation of the neutral stability curve (N), the Eckhaus (E) and 
Zig-Zag (ZZ) boundaries for the case of no flow. Within the shaded region the lamellar pat- 
tern is linearly stable. The effect of an imposed shear on a uniform lamellar structure can be 
qualitatively understood as the displacement of the state point along a line of constant e. This 
representation, however, fails to give the correct location of the stability boundaries given in the 
text. 
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FIG. 5. Stability diagram of a lamellar pattern under oscillatory shear of amplitude 7 = 0.2 
at e = 0.04. The solid lines show the neutral stability curve. The two inner lines marked with 
circles and squares bound the region of stability of the lamellar structure against long wavelength 
perturbations. 
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FIG. 6. Stability diagram of the lamellar pattern under oscillatory shear shear of amplitude 
7 = 0.4 at e = 0.04. The solid lines show the neutral stability curve. The two inner lines 
marked with circles and squares bound the region of stability of the lamellar structure against long 
wavelength perturbations. 
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FIG. 7. Results of the numerical integration of the model equation to show the instability of a 
lamellar pattern (shown in grey scale). The field ip shown in this figure has been transformed back 
to the laboratory frame of reference. Top, left to right: instability approximately of the zig-zag 
type, followed by kink band formation and reconnection leading to a different orientation. Bottom, 
instability approximately of the Eckhaus type leading to a wavelength increase without change in 
the orientation. 
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